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^T) ■ Cosmogenic isotopes ^Li and ^He produced in the detector are the most problematic 

^^ i background in the reactor neutrino experiments designed to determine precisely the 

^£ I neutrino mixing angle ^13. The average time interval of cosmic-ray muons in the 

\^ • detector is often on the order of the lifetimes of the ^Li and ^He isotopes. We have 

^^ , developed a method for determining this kind of background from the distribution 

^ [ of time since last muon for muon rate up to about 20 Hz when the background-to- 

I ■ signal ratio is small, on the order of a few percents. 
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1 Introduction 



Recently new experiments are proposed to determine the neutrino mixing an- 
gle ^13 using reactor anti-neutrinos with sensitivity in sin^26'i3 on the order of 
0.01 [1], about a factor of ten better than the current limit [2]. Comparing to 
the past reactor neutrino experiments [2,3,4], the systematic error of the pro- 
posed experiments will be greatly reduced by carrying a relative measurement 
with detectors located at two different baselines to cancel the reactor-related 
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errors, and other improvements to the detectors for reducing the detector sys- 
tematic errors. However, the cosmogenic ^Li and *He background will remain 
an irreducible systematic uncertainty in the measurement. 

The proposed reactor neutrino experiments detect the anti-neutrinos via the 
inverse /5-decay reaction with liquid scintillator. The inverse /3-decay process is 
the interaction of an electron anti-neutrino with a proton (hydrogen) resulting 
in the production of a positron and a neutron. A neutrino event is identified 
by the coincidence of the energy deposited by the positron (1 MeV to 8 MeV) 
followed by the energy released from neutron capture. Although the proposed 
experiments will be carried out underground, a nontrivial number of cosmic- 
ray muons can still pass through the hermetically shielded detectors. Some of 
these cosmic muons interact with the carbon nuclei in the liquid scintillator to 
produce unstable isotopes. Among them, ^Li and ^He have /9-neutron cascade 
decay modes, with branching fraction of 49.5% and (16 ± 1)%, respectively. 
The /9-decay of ^Li and ^He has a Q value of 13.6 MeV and 10.7 MeV, re- 
spectively, which overlaps with the positron signal of the neutrino events. The 
beta particle can thus mimic the prompt signal of a neutrino event while the 
neutron provides the delayed signal. Hence the cosmogenic ^Li and ^He iso- 
topes are indistinguishable from the neutrino events. Since these two isotopes 
are relatively long-lived, with a half-life of 0.178 s for ^Li [5] and 0.119 s for 
^He [6], they are difficult to reject with a muon veto. 

Production of ®He and ^Li has been measured with 190 GeV muons on a liquid- 
scintillator target at CERN [7]. Since their lifetimes are so close it is hard 
to determine their individual cross sections. The cross section ^(^Li + ^He) 
is measured to be (2.12 ± 0.35);ub, and the energy dependence of the cross 
section is suggested to be atot{E^) oc E^, with a = 0.73. In the KamLAND 
experiment, about 85% of the isotopes are found to be produced by inelastic 
scattering of energetic muons [4,8]. From the time distribution and the j3- 
energy spectrum of the /5-neutron events, the contribution of ^He relative to ^Li 
is less than 15% at the 90% confidence level. Furthermore, the ^He contribution 
could be determined by tagging the cascade ^He — i>*Li — i>®Be decay sequence 
[9]. Hence, for simplicity, we assume that the cosmogenic isotope background 
is solely ^Li in our investigation. 

While KamLAND can measure the ^Li background well from the time distri- 
bution of the /3-neutron event since last muon, it is not necessarily true for 
other experiments with less overburden. The average overburden of the Kam- 
LAND detector is 2700 m.w.e. (meter- water-equivalent), resulting in a muon 
flux of 0.0015 Hz/m^, or 0.2 Hz of cosmic-ray muons in the active volume of 
the detector. This implies the mean time interval of two adjacent muons is 5 s, 
much longer than the lifetime of ^Li. For most reactor experiments proposed 
to measure 6*13, the overburden ranges from about 300 to 1 200 m.w.e. [1]. 
The muon rate can be as high as 20 Hz, resulting a mean time interval of 



muons significantly shorter than the hfetime of ^Li. However, we find that it 
is still possible to determine the ^Li content since the anti-neutrino rate and 
the ^Li rate are significantly lower than the muon rate. The probability of 
having two /5-neutron events (either neutrino or ^Li candidates here) in the 
time window used in the analysis is negligible. Thus the muons immediately 
following neutrino candidates can simply be ignored. In this paper, we show 
that it is possible to measure the ^Li and ^He background in-situ even when 
the cosmic-ray muon rate, comparing to the isotope lifetime, is relatively high. 



2 Method of Least Squares 



Suppose the muons are randomly distributed in time. The time interval t of 
two adjacent muons follows a normalized exponential function 
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where the constant T is the mean time interval between two adjacent muons, 
T = l/Rn, with Rn being the muon rate. The neutrino events are also ran- 
domly distributed and independent of the muons. The time interval between 
a neutrino event and the last muon preceding it obeys the same exponential 
distribution, but with a different constant T' = l/(-R^ + Ru), where R^, is 
the neutrino rate. The neutrino rate is normally several orders lower than the 
muon rate. The difference in T' and T can thus be ignored. Therefore, we have 
fuit) = f]i{t). Here, the variable t is the time of the neutrino event since the 
last muon. 

The ^Li events are correlated with the muons. If the muon rate is high, a 
^Li event might be the product of any of the preceding muons. When summed 
over the contributions of all preceding muons, the probability density function 
(p.d.f.) of the ^Li events in terms of the time since last muon t is 
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where r = 0.257 s is the lifetime of ^Li, and ti is the time interval of two 
adjacent preceding muons. The first term corresponds to the contribution from 
the last muon. The second term is that from the second-to-the-last muon, and 



so on. Because the time since last muon is t, the time interval between the 
two muons immediately before and after a ^Li event must be greater than t. 
This condition has a probability of exp(— t/T). Therefore the p.d.f. of the time 
since last muon t for the ^Li events becomes 

A.(t) = ^exp(-t/A), \ = l + J^- (3) 



In a reactor neutrino experiment, the observed distribution of the time since 
last muon t for all /9-neutron events is a combination of the distribution of the 
genuine neutrino events and that of the ^Li events. As a result, the observed 
distribution can be described by 

fit) = B-j exp(-t/A) +S-^ exp(-t/T) , (4) 

where B and S are number of ^Li and neutrino events, respectively. The 
background-to-signal ratio is simply B/S. The amount of the ^Li background 
can then be measured by minimizing the x^ function which is a sum of the 
squares of the residuals of the observed time distribution divided in bins from 
the expectation given by Eq. 4. 



3 Method of Maximum Likelihood 



In the least-squares method, contributions of all preceding muons are included 
in the p.d.f. and expressed as a function of the time since last muon. Similarly, 
an unbinned maximum-likelihood fit can be performed with the same p.d.f.. 
The log-likelihood function (MLa) is thus defined as 



logL = ^log l^ie-*'/^ + (1 - b)^e-'^/^' 



(5) 



where tj is the time of the i-th /3-neutron event since last muon. Here we 
assume that each /3-neutron event has a probability of b being a ^Li event and 
1 — 6 to be a neutrino candidate. To precisely determine ^13, it is desirable to 
have the background-to-signal ratio as small as possible. In this case, b can be 
taken as such a ratio. 

Alternatively, the maximum-likelihood method can exploit the contributions 
to a /9-neutron event from each preceding muons, instead of averaging them 
into the time distribution since last muon. When the muon rate is not too high 



but contributions from the preceding muons other than the last one cannot 
be ignored, we can construct a log-hkehhood function (ML) as 
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where i sums over all /9-neutron events in the sample and j sums over all 
preceding muons of the i-th /5-neutron event. Here tn is the time of the i-th 
/9-neutron event since last muon, and tij is its time since the j-th preceding 
muon. As in the least-squares method, the term exp(— tji/T) is included in 
the probability density function of the ^Li to account for the condition that 
there is no muon between the /^-neutron event and the last muon. In practice, 
only muons in a finite time window are included. This truncation will lead to 
an error in the probability on the order of 10~^ for a two-second window. 

The statistical error of the parameter in the fit can be estimated for the 
maximum-likelihood method when the sample size is very large. In this case, 
the variance V of the parameter 9 is given by [10] 
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where / is the density function, 9 is the expectation of the parameter ^, X is 
the sample space, and E denotes the mathematical expectation in the sample 
space. For the MLa given in Eq. 5, the variance is 
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If 6 ^ r/(r + T), which is the case for most reactor neutrino experiments, the 
above integral can be simplified. The statistical error of b is approximated to 
be 
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where N is the total event number of /5-neutron events. It is interesting to note 
that the error is independent of b. When the muon rate is high, as expected, 
the error gets worse. 



4 Monte Carlo Simulation 



To validate the formulation, we have simulated data samples of a reactor 
neutrino experiment, with neutrino, muon, and ^Li events taken into account. 
For each sample, two independent data sets are first generated, one contains 
neutrino events, and the other contains muons and ^Li events. The neutrino 
events are generated randomly in time. For a given muon rate i?^, the number 
of muons is determined. Muons are then generated randomly in time. The 
probability of producing the ^Li isotope by each muon is calculated with a 
Poisson distribution based on the production cross section measured at 190 
GeV. If a ^Li is created, its decay time is generated according to an exponential 
function using the ^Li lifetime. Neither the energies of the ^Li and neutrinos 
nor the muon-energy dependence of the yield of ^Li is taken into account here. 
The two independent data sets are then combined and sorted by time to form 
the /3-neutron sample for analysis. 

The distribution of the time since last muon for all the /?-neutron events in a 
Monte Carlo sample with an input muon rate of 4 Hz and a value of 0.01 for 
B/S is shown in Fig. 1. The red curve is the best fit of Eq. 4 to the distribution, 
whereas the blue curve is a fit with B fixed to and S as a fit parameter. Since 
the amount of ^Li background is small, the two curves are almost identical. 
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Fig. 1. Distribution of time since last muon for Monte Carlo simulated /3- neutron 
events (a combination of genuine neutrino and ^Li events). The input muon rate is 
4 Hz and the B/S ratio is set at 0.01. The red line is the best fit of Eq. 4 to the 
distribution. The blue line is a fit with B = and S as a fit parameter. 

The results obtained with the least-squares and maximum-likelihood (ML) 
methods are shown in Fig. 2. For each given muon rate, the background-to- 



signal ratio, B/S, is fixed at 1% since the statistical error is independent of it 
(see Eq. 9). In addition, 400 samples are generated and analyzed for each given 
muon rate to investigate the potential bias and precision of the fit. It is clear 
from Fig. 2 that both the least-squares and maximum-likelihood methods can 
determine the background-to-signal ratio reliably up to a muon rate of about 
20 Hz. 
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Fig. 2. Background-to-Signal ratios (B/S) obtained with the least-squares and 
maximum-likelihood methods as a function of the muon rate. The input back- 
ground-signal ratio is fixed at 1%. For each muon rate, 2.5 x 10"' neutrino 
events are generated, corresponding to a 0.2% statistical error. The results from 
the least-squares fit are plotted to the right of those obtained with the maxi- 
mum-likelihood fit for clarity. 

Along with the analytical estimation of Eq. 9, the relative resolution of B/S 
obtained with the least-squares and maximum-likelihood fits as a function of 
the muon rate are shown in Fig. 3. Here the relative resolution is defined as the 
statistical error (standard deviation) obtained from the fit normalized to the 
input value of 1% for the B/S ratio. The results from the two fitting methods 
are in excellent agreement with the analytical calculation. 



5 Conclusions 



We have shown that, as long as the occurrence of (anti-)neutrino events and ^Li 
decays is low comparing with the ^Li lifetime, it is possible to measure in-situ 
the cosmogenic ^Li background for a reactor neutrino experiment with rela- 
tively high muon rate. The muons immediately following a /3-neutron event 
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Fig. 3. Relative resolution of B/S as a function of the muon rate, The solid line is 
the analytical result of Eq. 9. The relative resolution is defined as the statistical 
error (standard deviation) obtained from the fit normalized to the input value of 
1% for the B/S ratio. 

can be safely ignored because they have very small probability to be asso- 
ciated with another /3- neutron event. We have obtained probability density 
functions that can be utilized for determining the amount of ^Li background 
using the least-squares or maximum- likelihood fit, with the latter method hav- 
ing a slightly better precision. When the background-to-signal ratio b is small 
such that b <C t/{t + T), the precision of the in-situ measurement can be 
estimated with a simple formula in Eq. 9, which depends on the muon rate 
and the number of /3-neutron events in the sample. 
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